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Abstract. Under a convexity assumption on the boundary we solve a local 
inverse problem, namely we show that the geodesic X-ray transform can be 
inverted locally in a stable manner; one even has a reconstruction formula. We 
also show that under an assumption on the existence of a global foliation by 
strictly convex hypersurfaces the geodesic X-ray transform is globally injective. 
In addition we prove stability estimates and propose a layer stripping type 
algorithm for reconstruction. 



1. Introduction 

Let X be a strictly convex domain in a Riemannian manifold (X, g) of dimension 
> 3. In this paper we consider the local inverse problem for the geodesic X-ray 
transform. That is, for an open set O C X, we call geodesic segments 7 of g which 
are contained in O with endpoints at dX O-local geodesies; we denote the set of 
these by M.o- Thus, Aio is an open subset of the smooth manifold of all geodesies, 
M. . We then define the local geodesic transform of a function / defined on X as the 
collection (If)(j) of integrals of / along geodesies 7 £ M-O: i- e - as the restriction 
of the X-ray transform to A4q- 

In order to state our main theorem in concrete terms, it is useful to introduce 
some notation. Let p £ C°°(X) be a defining function of dX, considered a function 
on X (so p > in X, < on X \ X, vanishes non-degenerately at dX). Our main 
theorem is an invertibility result for the local geodesic transform on neighborhoods 
of p in X of the form {x > — c}, c > 0, where J is a function with x{p) = 0, 
dx(p) = —dp{p), see Figure [I] below. 

Theorem. For each p £ dX, there exists a function x £ C°°(X) vanishing at p 
and with dx(p) = —dp(p) such that for c > sufficiently small, and with O p — {x > 
— c} H X, the local geodesic transform is injective on H s (O p ), s > 0. 

Further, let H s (Mo p ) denote the restriction of elements of H S (M) to A4o p , 
and for f > let 

H s F {O p ) = e F '^H s = {f£ H? oc (O p ) : e -^+ c )/ £ H s {O p )}. 

Then for s > there exists C > such that for all f £ H?r(O p ), 

ll/Hif'-^Op) ^ C 'll / /l>lo p ll^(At 0p )- 
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Remark. Here the constant C is uniform in c for small c, and indeed if we consider 
the regions {p > Pq} D {x > — c} with \p \ and |c| sufficiently small and such that 
this intersection is non-empty, the estimate is uniform in both c and pq. 

Further, the estimate is also stable under sufficiently small perturbations of the 
metric g, i.e. the constant is uniform. (Notice that the hypotheses of the theorem 
are satisfied for small perturbations of g\) 

We remark that for this result we only need to assume convexity near the point 
p. This local result is new even in the case that the metric is conformal to the 
Euclidean metric. We also point out that we also get a reconstruction method in 
the form of a Neumann series. See Section [5] for more details. 

While this large weight e F ^ x+c ^ means that the control over / in terms of // 
is weak at x — — c, the control is uniform in compact subsets of O p : these weights 
are bounded below on O p by a positive constant, and bounded above on compact 
subsets of O p (in particular at parts of dX). Here F > can be taken small, but 
not vanishing. Further, x, whose existence is guaranteed by the theorem, is such 
that x = — c is concave from the side of O p . 

As an application, we consider domains with compact closure X equipped with 
a function p : X —} [0, oo) whose level sets E t = p _1 (i), t < T, are strictly convex 
(viewed from / o~ 1 ((i, oo)) (and dp is non-zero on these level sets), with So = dX 
and X\U te [o i T)^t — P~ {[T, oo)) either having measure or having empty interior. 
(Note in particular that p is a boundary defining function.) 

Corollary. For X and p as above, if X \ U t6 [ 0j T)^t nas measure, the global 
geodesic transform is injective on L 2 (X), while if it has empty interior, the global 
geodesic transform is injective on H S (X), s > n/2. 

This corollary is an immediate consequence of our main theorem. Indeed, if 
If = and / E H s , s > n/2, f ^ 0, then supp/ has non-empty interior since / 
is continuous by the Sobolev embedding, while if / <E L 2 , f ^ 0, then supp / has 
non-zero measure. On the other hand, let r = inf supp f p; if t < T we are done, for 
then supp / Cl \ U te [o.T)St. Thus, suppose r > T, so / = on S t for t < r, but 
there exists q € S T n supp / (since supp / is closed and X is compact) . Now we 
use the main theorem on p _1 (r, oo) to conclude that a neighborhood of q is disjoint 
from supp / to obtain a contradiction. 

In fact, in this global setting we can even take x = —p, and the uniformity 
of the constants in terms of c and po, as stated in the remark after the main 
theorem directly yields that if t < T then there exists S = 6t > such that if c, po € 
(t — S t ,t+6t) then a stability estimate holds (with a reconstruction method!) for the 
region p^ 1 ([po, c)). Now in general, for T' < T , one can take a finite open cover of 
[0, T'] by such intervals (t'pt"), j = 1, . . . , k (with, possibly after some reindexing 
and dropping some intervals, t[ < 0, t'^ > T' , t" 6 (lfj +1 ,t" +1 )), and proceed 
inductively to recover / on U te [o,T']£t from its X-ray transform, starting with the 
outermost region. More precisely, first, using the theorem, one can recover the 
restriction of / to p _1 ((— oo, t'{)). Then one turns to the next interval, (t^t^), and 
notes there is a reconstruction method for the restriction to /0 _1 ((*2; ^2)) 01 functions 
fi supported in p~ x {(t' 2 , +00)) (no support condition needed at the other end, t'2). 
One applies this to f% = fcf, where (j>2 identically 1 near p _1 ([t", +00)), supported 
in p^ 1 {{t' 2 , +<x>)); since / = (1 — 02)/+<fo/) an d one has already recovered (1 — 02)/, 
one also knows the X-ray transform of 2 /, and thus the theorem is applicable. 
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One then proceeds inductively, covering p~ ([0,T']) in k steps. This gives a global 
stability estimate, and indeed a reconstruction method doing a reconstruction layer 
by layer; that is, we have (in principle) developed a layer stripping algorithm for 
this problem. 

The geodesic ray transform is closely related to the boundary rigidity prob- 
lem of determining a metric on a compact Riemannian manifold from its boundary 
distance function. See [T71I3] for recent reviews. The case considered here is the lin- 
earization of the boundary rigidity problem in a fixed conformal class. The standard 
X-ray transform, where one integrates a function along straight lines, corresponds 
to the case of the Euclidean metric and is the basis of medical imaging techniques 
such as CT and PET. The case of integration along more general geodesies arises in 
geophysical imaging in determining the inner structure of the Earth since the speed 
of elastic waves generally increases with depth, thus curving the rays back to the 
Earth surface. It also arises in ultrasound imaging, where the Riemannian metric 
models the anisotropic index of refraction. Uniqueness and stability was shown by 
Mukhometov on simple surfaces, and also for more general families of curves in 
two dimensions. The case of geodesies was generalized also for simple manifolds to 
higher dimensions in [13] . pQ. In dimension n > 3, the paper [3] proves injec- 
tivity and stability for the X-ray transform integrating over quite a general class of 
analytic curves with analytic weights, assuming an additional microfocal condition 
that includes the case of real-analytic metrics for a class of non-simple manifolds. 
Reconstruction procedures or inversion formulas have not been proven except in 
a few cases for instance for a class of symmetric spaces, see [5], and real-analytic 
curves [3J . Our results generalize support type theorems to the smooth case for the 
geodesic X-ray transform given in [5] for simple real-analytic metrics. 

The global geometric condition that we are imposing is a natural analog of the 
condition -ip(r/c(r)) > proposed by Herglotz [6] and Wiechert and Zoeppritz [21] 
for an isotropic radial sound speed c(r). In this case the geodesic spheres are strictly 
convex. It is also satisfied for negatively curved manifolds. But this condition allows 
in principle for conjugate points of the metric. In |19) one can find a microfocal 
study of the geodesic X-ray transform with fold caustics. A similar condition of 
foliating by convex hypersurfaces was used in |20j to satisfy the pseudoconvcxity 
condition needed for Carleman estimates. 

We also remark that our approach is a completely new one to uniqueness for the 
global problem for the geodesic ray transform. The only method up to now, except 
in the real-analytic category [17], has been the use of energy type equalities one 
introduced by Mukhometov [11 and developed by several authors which are now 
called "Pestov identities" . 

The main theorem is proved by considering an operator A which is essentially a 
'microfocal normal operator' for the geodesic ray transform. Let p be a boundary 
defining function of X, i.e. p > in X, p = at dX, and dp ^ at dX; we 
assume that in fact p is defined on the ambient space X as above. First we choose 
an initial neighborhood U of p in X and a function x defined on it with x(p) = 0, 
dx(p) — —dp(p), dx ^ on U with convex level sets from the side of the sublevel sets 
and such that O c = {x > — c}P\{p > 0} satisfies O c C U is compact. Such a x exists 
as can be seen by slightly modifying — p, making the level sets slightly less convex. 
We define an operator L which integrates If over a subset of M.o a with a C°° cutoff, 
and consider A = L o I. We consider this operator as a map between appropriate 
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function spaces on O c . It turns out that with the subset of geodesies we choose, 
the exponential conjugate Af of A is a pseudodifferential operator in Melrose's 
scattering calculus (10j . (The exponential conjugate corresponds to working with 
exponentially weighted spaces for A.) We show that Af is a Fredholm operator, 
and indeed that it is invertible for c near 0. 

X 



x 





Jp ^^^^ 










= 


x = -^c\ 

,p X 



Figure 1. The functions p and x when the background is flat 
space X. The intersection of p > and x c > (where x c = 
x + c, so this is the region x > — c) is the lens shaped region O p . 
Note that, as viewed from the superlevel sets, thus from O p , x has 
concave level sets. At the point z, L integrates over geodesies in 
the indicated small angle. As z moves to the artificial boundary 
x c = 0, the angle of this cone shrinks like Cx c so that in the limit 
the geodesies taken into account become tangent to x c = 0. 



Before giving more details, recall that Stefanov and Uhlmann jTSj have shown 
that under a microlocal condition on the geodesies, one can recover the singularities 
of functions from their X-ray transform, and indeed from a partial X-ray transform 
(where only some geodesies are included in the X-ray family A4'). (In fact, they 
also showed analogous statements for the transforms on tensors.) Roughly speaking 
what one needs is that given a covector v — (z, C), one needs to have a geodesic 
in M! normal to £ at z such that in a neighborhood of v a simplicity condition 
is satisfied. Indeed, under these assumptions, a microlocal version of the normal 
operator, (QI)*(QI), where Q microlocalizes to M! roughly speaking, is an elliptic 
pseudodifferential operator. Now, in dimension > 3, if the boundary OX is convex, 
one can use geodesies which are almost tangent to dX to give a family M! which 
satisfies the above conditions for v with z near dX. While this gives a recovery 
of singularities for the local problem we are considering, it yields no invertibility 
or reconstruction. Indeed for the latter we would like to have (QI)*(QI) to be an 
invertible operator on a space of functions on O c ; in particular, as one approaches 
x = — c one would need to only allow integrals over geodesies in a narrow cone, 
becoming tangent to x — — c, which takes one outside the framework of standard 
pseudodifferential operators. 
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To remedy this, we introduce the artificial boundary x = — c, and work with 
pseudodifferential operators in x c = x + c > which degenerate at x c — 0. Sup- 
pressing the c dependence of x, the particular degeneration we end up with is 
Melrose's scattering calculus as already mentioned. This is defined on manifolds 
with boundary, with boundary defining function x, and is based on degenerate vec- 
tor fields x 2 d x and xd Vj , where the (x, y\, . . . , y n -i) are local coordinates. This 
has the effect of pushing x = 'to infinity' (these vector fields are complete un- 
der the exponential map). Thus, ultimately, our approach is based on working in 
a framework with an artificial boundary which is effectively 'at infinity', and we 
work with function spaces allowing exponential growth at this boundary. Thus the 
control at x = will be quite weak in a sense, though one has the standard control 
when x is bounded away from 0. Since x = is just an artificial boundary, this is 
a satisfactory situation. 

In fact, for most of the paper we work in a much more general setting. We 
consider a family of curves 7^ : / — > X parameterized by v = (z, () e SX (the 
sphere bundle of X realized as a subbundle of TX, e.g. via a Ricmannian metric) 
with 7^(0) = v and we assume that if v is tangent to a level set of x in O c , i.e. if 
^(x o 7„)| t =o = 0, then ^(x o 7 i ,)| t=0 > C > 0. By possibly shrinking U, we may 
always assume this in our setting; the lower bound on the second derivative is a 
concavity statement for the level sets of x from the side of the superlevel sets. Let 
x = x c = x + c as above. Thus, x is a boundary defining function for {x > — c}; 
for the time being we regard c as fixed. A consequence of our uniform concavity 
statement is that, with A = ^(x o 7„)| t= o, if C\ > is sufficiently small and 
|A| < Ciyfx, then 7„ remains in x > 0. Rather than using this range of A, wc 
instead use the stronger bound |A| < C2X, and define A to be an average: 

Af(z)=x- 1 J If{ lv ) X {\/x)d l x{v), 

where jj, is a non-degenerate smooth measure on SX, and \ nas compact support. 
We show that for FeE, 

Af = x- 1 e- F > x Ae F ' x e > 0}), 

where stands for the scattering calculus of Melrose, and is elliptic in the sense 
that the standard principal symbol is such near the boundary (up to the boundary, 
x = 0). However, even when this holds globally on a compact space, this ellipticity is 
not sufficient for Fredholm properties (between Sobolev spaces of order shifted by 1), 
or the corresponding estimates, due to the boundary x = 0. In general, scattering 
pseudodifferential operators also have a principal symbol at the boundary, which 
is a (typically non-homogencous) function on a cotangent bundle; this needs to be 
invertible (non-zero) globally to imply Fredholm properties. Similarly, estimates 
implying the finite dimensionality of localized (in O) non-trivial nullspace as well 
as stability estimates, follow if this principal symbol is also invertible on O. (Note 
that here localization does allow the support in {x > 0} to include points at x — 0!) 
We thus show that in the case when 4^ (oio7„)|j_o is a quadratic form in £ subject to 
A = 0, which is the case with geodesies, for suitable choices of x, namely essentially 
cutoff Gaussians, this principal symbol is invertible when the weight F satisfies 
F > 0. This implies that Af is Fredholm on this space, i.e. A itself is Fredholm on 
exponentially weighted spaces, where exponential growth is allowed at x — 0. We 
now recall that x — x c depends on c, with all estimates uniform for c remaining in 
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a compact set, and the argument is finished by showing that for c > sufficiently 
small one not only has Fredholm properties but also invertibility, essentially as the 
Schwartz kernel has small support. 

We note that the geodesic nature of the curves was only used in the crucial 
step of showing that the principal symbol at the boundary is invertible. While our 
argument relied on properties of the geodesies to analyze this symbol, it may well 
be possible to analyze it in general and prove the result for more general families 
of curves. We remark that J. Boman has given in [5] counterexamples for local 
uniqueness for the X-ray transform that integrates along lines with a dense family 
of smooth weights so that we expect some restrictions on the family of curves. 

2. Scattering calculus 

Melrose's algebra of scattering pscudodiffcrcntial operators ^™^(M) on a com- 
pact manifold with boundary M, see |10j . can be thought of either via reducing to 
a model on K n (via appropriate charts on M, the interior of M ), or via a geometric 
definition. Both are of use in the current paper; the K n version makes the simplicity 
of this algebra transparent, while the geometric definition emphasizes that infin- 
ity in the M™-picture is not really 'remote', and indeed in our setting the artificial 
boundary x = — c plays dM, i.e. infinity is at a decidedly finite place (moving it to 
infinity is what is artificial). 

First we start with the K™ picture, which is straightforward. Indeed, the scat- 
tering algebra in this setting is a special case of Hormander's Weyl calculus [7J 
Section 18.5], which in this particular case has also been studied by Parenti [T3] 
and Shubin [15]. That is, scattering symbols of order (to, I) are defined to be 
functions on M™ x K™ satisfying 

\D?D*a(z,0\ <C af3 (z) l -M(O m - m , 

i.e. they are 'product type' symbols in z and £. Note that our order convention 
for the second order Z, indicating growth/decay in z, is the opposite of that of 
Melrose [TU] (i.e. our I is — / in HO]); we make this deviation so that the symbol 
class increases both with to and I, i.e. so that the two indices play a parallel role. 
Their set is denoted by S m ' l (W n ,W n ) or simply S m > 1 . One then defines ^™> l (R n ) 
to consist of, say, left quantizations of such symbols, i.e. of operators of the form 

(2.1) A„(,) = ( 2 .)-»/ e «.-.'«„(,, fW y),^, 

understood as an oscillatory integral. Right quantizations could be used equally 
well, i.e. one gets the same class of operators if a £ S m ' but one substitutes a(z' , £) 
into the oscillatory integral in place of a(z, (). Note that for I < 0, ^™ ,z (K n ) is 
a subspace of Hormander's uniform algebra ^"(M"), i.e. where the above esti- 
mates hold without the factor (z)'~l Q l, and the general weight barely affects the 
standard arguments with pseudodifferential operators. The space ^^(R") is a 
filtered *-algcbra under composition of operators and taking adjoints (relative to 
the Euclidean metric), i.e. 

A £ * s " c M (ir), B £ *™'' r (M") => AB £ *™+ m '''+ r (M™) 

and 

A £ *^(R n ) => A* £ ^''(R n ). 
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Further, we define the principal symbol of A to be the equivalence class of the am- 



plitude a in ([2T) in S m ' 1 / , which thus captures A modulo f "-L'-^r), 
i.e. one order lower operators both in terms of the differential order and growth 
at infinity. With this definition, the principal symbol of AB is the product of 
the principal symbols of A and B, while that of A* is the complex conjugate of 
the principal symbol of A. In particular, if A is elliptic, i.e. its principal sym- 
bol is invertible in the sense that there is b € S~ m ' such that ab — 1 € 
(which is independent of the choice of representative for the principal symbol), 
then the standard parametrix construction produces B € \E'~ m ' _i (R n ) such that 
AB - Id e t^ 00 ^ 00 ^ 11 ). Operators R in * E T C °° -°°(R n ) have a Schwartz function 
on R 2n for their Schwartz kernel; this is just the inverse Fourier transform of their 
amplitude r in the ( variable evaluated at z — z' (where R 2 ™ — R™ x R™, , with z the 
left and z' the right variable) . In particular, such operators are compact between all 
polynomially weighted Sobolev spaces H s < r = (z)- r H s (W l ). Further, A G *™' ; (R n ) 
is bounded H s ' r —} H s ~ m ' r ~ l , and if A is elliptic then the parametrix construction 
and the compactness we observed shows that A is Fredholm - it has closed range, 
finite dimensional kernel and cokernel, and corresponding estimates, 

IMIff- < C(\\A U \\ H s-m,r-l + \\FU\\ H -N,-N), 

where F can be taken a finite rank element of 4'~ 00 '~ 00 (IR n ), and N can be taken 
arbitrary. 

In order to relate \I/ sc (R n ) to the geometric setting, and also in order to explain 
its classical subalgebra, it is useful to compactify R™. Concretely, we compactify 
1" to a closed ball 1™ by adding the sphere at infinity S™" 1 . Thus, R" \ {0} can be 
identified with (0,oo),. x Sg 1_1 via 'polar coordinates', (r, 0) i-> r8; letting x = r _1 
we have 'reciprocal polar coordinates', (0, oo) x x Sg -1 which allow us to glue a 
sphere to x — (corresponding to r — oo) by extending the range of x to [0, oo). 
(Thus, formally, M." is the disjoint union of R" with [0, oo) x §™ _1 modulo the 
identification of R" \ {0} with (0,oo) a x Sg -1 .) Notice that x — r 1 is a boundary 
defining function near <91R™; modifying it near gives a global boundary defining 
function p. It is straightforward to check that Schwartz functions on R ra are exactly 
the restrictions to R" of C°° functions on R n which vanish with all derivatives at 
9R™. Further, writing z as the variable on R", the linear vector fields Zjd Zk on R n 
lift (automatically uniquely, as R™ is the interior of R") to smooth vector fields on 
R n which are tangent to the boundary, and indeed all smooth vector fields tangent 
to the boundary are, away from the origin, linear combinations of these lifts with 
coefficients that are smooth on R™. Since being a symbol on R™, i.e. satisfying 
estimates |-D™a(z)| < C Q (z)'~' Q ', is equivalent (away from the origin, near which 
one has smoothness) to satisfying stable estimates under linear vector fields, i.e. 
that |Vi . . . Vfca| < C(z) 1 for all k and linear vector fields Vj (with C depending on 
these), it follows that the lift of a symbol is a conormal function, i.e. a function 
that satisfies pV\ . . .Vua £ L°° whenever Vj are vector fields tangent to cM™, 
and conversely, every conormal function is the lift of a symbol. Correspondingly 
/9~'C°°(R n ) C ^'(R™); these are the 'classical' or 'one-step' symbols; the Taylor 
series of a C°° function at the boundary gives rise to the expansion (with x = p 
near x = 0) 



E 
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understood as an asymptotic sum. 

One can now compactify each factor of R™ x R™ to define the compactified space 
of scattering symbols R™ x R"; we write pg for the boundary defining function in 
the first factor ('position', z) and p^ for that in the second factor ('momentum', 
Q). The same considerations as above show that a scattering symbol on R™ x R" 
of order (m,l) corresponds to a conormal function on R™ x R n , i.e. one satisfy- 
ing p^PgVi . . . Vka e L°° whenever Vj are vector fields tangent to both bound- 
ary hypersurfaces of R" x R™. Classical symbols, as before, then arc elements of 
Pg V^ m C°°(R" x R"), i.e. functions of the form a = p d l p^ m a, a e C°° (R™ x R" ) . 
Note that for a classical symbol, its equivalence class in S m ' 1 / S™ 1 ^ 1 ' 1 ^ 1 can be rep- 
resented by Pq V^o" 1 times the function ciq = a|,g(i^xl^) on 9(R n x R") in the sense 
that any smooth extension a' of this function to R" x R" produces an element of 
the equivalence class of a. Ellipticity then simply means the non- vanishing of this 
function a^. Note also that this principal symbol can be thought of as consisting 
of two parts, namely the standard principal symbol, at p^ = 0, and the 'boundary 
principal symbol' at p d = 0. We also write #™''(R") = * SC (R"). 

If M is a manifold with boundary with interior M, we can now define ^"''(M), 
much as the standard pseudodifferential algebra is defined on manifolds by locally 
identifying the manifold with R n and imposing that on such charts U x U the 
Schwartz kernel of the operator is that of a pseudodifferential operator on R", and 
allowing additional globally smooth terms in the Schwartz kernel. In our case, the 
analogous construction is locally identifying M with R", and imposing that on such 
charts U x U the Schwartz kernel of the operator is that of an element of ^"''(R"), 
and allowing additional globally Schwartz (i.e. rapidly decaying with all derivatives, 
smooth) terms in the Schwartz kernel. As in the standard manifold case, all the 
basic properties of the algebra generalize (one needs to impose some proper support 
conditions in the absence of compactness) . Concretely, the weighted Sobolev spaces 
H*f(M) are also defined by local identification with and then A e *™-'(R") 
implies that A is bounded from H^ r (M) to H^~ m ' r ~ l (M). 

It is also of some use to work out the behavior of the Schwartz kernel of elements 
of ^™' l (M) on M x M . In view of the previous definition, this reduces to a calcu- 
lation for f™ , '(IR' 1 ) (modulo Schwartz terms which we ignore as they give elements 
of \E , ~ 00, ~ 00 (M)). Thus, using local coordinates y on § n_1 , and corresponding 
coordinates (x,y,x',y') on M x M, one checks that in the coordinates 

x, y, x = — — , y = , 

X z X 

valid for x > 0, so the diagonal is X = 0, Y = when x > 0, the Schwartz kernel of 
an element of "3/™'' is of the form x~ l K, where K is smooth in (x, y) down to x = 
with values in conormal distributions on R^ Y , conormal to {X = 0, Y = 0}, which 
are Schwartz at infinity (i.e. decay rapidly at infinity with all derivatives). Further, 
the boundary principal symbol is simply x~ l times the Fourier transform in (X, Y) 
of K\ x= o (a restriction which makes sense in view of the stated smoothness). In 
particular, when / — 0, we need to check that 



/ 



e- i£x - ir >- Y K{0, y, X, Y) dX dY 



is a non-zero function of (y,£,i]), with a lower bound C((^,ri)), C > 0, for its 
absolute value (which means we also need a uniform bound at infinity in addition 
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to the invertibility) . Checking this will be the main step of the arguments presented 
in the next section. 

We mention here that vector fields in V sc (Af) = xVb(M), where Vb(M) is the 
set of all smooth vector fields tangent to dM, are in ^lf(M), and indeed the 
Sobolev spaces of positive integer differential orders s are equivalently defined by 
u £ H*f (M) if and only if x~ r Vi . . . V k u £ L 2 C (M) for all k < s (including k = 0) 
and Vj £ V SC (M); here L 2 C (M) is the L 2 space given by identification by K n , i.e. 
the measure (or density) is, up to a non-degenerate positive multiple, r™ -1 drdy = 
x~ n ~ x dxdy. (Densities like this may be called scattering densities.) 

We now briefly relate the standard Sobolev spaces H S (M) to H^ r (M) for s > 0. 
First, for s = 0, the above description gives #s°<r ( " +1)/2 (M) = L 2 (M) (in the 
sense of equivalent norms). Next, using that V[ . ■ ■ V fc 'it £ L 2 (M) for k < s and 
Vj smooth vector fields on M (which is equivalent to u £ H S (M)) implies that 

V x ...V k u £ if s ° c '" (n+1)/2 (M) for k < s and V 3 £ V SC (M) (since all elements of 
Vsc(M) are smooth vector fields), i.e. that u £ H%£~^ n+1 ^ 2 (M), so 

n 4- 1 

(2.2) H S (M) C H% r {M), r< 

with continuous inclusion map. For the converse direction, we note that if V is a 
smooth vector field, then x 2 V £ V SC (M). Thus, V x . . . V k u £ H°i 2s ~ (n+1)/2 (M) for 
k < s and Vj £ V SC (M), so x~ 2s Vi ...V k u£ fl s ° c " ( " +1)/2 (M), so x" 2 ^ . . . X - 2 V k u £ 
Hsi~ (n+1)/2 (M), implies that V{ . . . V^u £ L 2 (M) for k < s and Vj smooth vector 
fields. Thus, 

(2.3) H^(M) C H S (M), r> ^ + 2s, 

with continuous inclusion map. There are similar inclusions between negative order 
spaces. For instance, as H~ S {M) — (Hq(M))*, s > 0, via identification by the L 2 



pairing, and as Hq(M) is a closed subspace of H S (M), the inclusion (2.2 1 gives the 
continuous inclusion map on the dual spaces 

(2.4) H-*>- r (M) c H-°(M), -r > 

Finally we discuss what happens when ellipticity holds only locally. Thus, sup- 
pose O is an open subset of M on which A £ t&™> l (M) is elliptic, and suppose 
that K C O is a compact subset. Let <p be supported in O, identically 1 on K; 
let O' be a neighborhood of supp</> with closure compactly contained in O. By the 
ellipticity assumption, there is a local parametrix G £ 1 I r ~ m— ' (M) for A such that 
GA = Id+E, E £ *°^°(M), but over O' the better conclusion that E is, locally, 
in * sc 00 '" 00 , holds, so fiEfi £ -°°(M). Thus, (j)E(j) is compact on any poly- 

nomially weighted Sobolev space, so in particular there is a finite rank operator 
F £ * SC 00 ' _00 (M) supported in O x O such that Id+cpEcp - F is invertible. Now 
suppose that v is supported in K, so cf>v = v. Then <fiGA(f> = <p 2 + tf>E<j) shows that 

(ld+cj)E(t>)v = (f>GAv, 

so 

v = (Id +<f>E(j) - F)~ 1 4>GAv - (Id +(j)E(j) - F)~ 1 Fv. 

In particular, if Av = then v is in a finite dimensional space, namely the range of 
(Id +<j)E4>— F)~ 1 F, and if one chooses a complementary subspace V of Ker AC\ { 



w 
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suppw C K} in a weighted Sobolev space H^ r (M), then there is a constant C > 
such that for v £ V n {w : supp to C If}, 

IMIff s r(M) ^ C 'll j4w llff B » c - ro ''- , (M)' 

i.e. a stability estimate holds. 

Now suppose that one has a family of operators, At, t £ [0,T], depending con- 
tinuously on t in &™> l (M), with each element of the family being elliptic on O 
(and thus there is a uniform constant in the estimates over compact subsets of 
O). Suppose also that we have a continuous function / on [0,T] with /(0) = 0, 
a compact subset K of O, and a family of open sets O t , t > in K , with the 
boundary defining function satisfying x < /(i) on O t , and we are interested in 
distributions v supported in Ot- In view of the uniform elliptic estimates, choosing 
O' a neighborhood of K with closure compactly contained in O, we then have fam- 
ilies of operators G t and E t , depending continuously on t £ [0, T], with values in 
y- m '~ l (M), resp. y° s i°(M), such that on O 1 , E t is uniformly in -°°(M). Thus, 
the Schwartz kernel K t of satisfies that for any N, x~ N (x')~ N K t is bounded 
(with values in scattering densities in the right, i.e. primed, factor), i.e. locally is 
of the form n t (y with \n t (x, y, x' , y')\ < CnX N (x') N . (Notice that the fact 
that we used 'scattering' densities is thus of little relevance; any polynomial factor 
such as (x')~™ _1 , can make no difference.) If <j> t £ C^°(M) is now supported in O t 
and takes values in [0,1], then (f> t E t 4>t has kernel <j>t{x, y)(j)t{x' , y')nt (Jnn+i , with 
\cpt(x,y)<i)t(x',y')Kt\ < C' N f(t) 2N x n+1 (x') n+1 for all N, and thus by Schur's lemma 
is bounded on L 2 C (M) with norm < C'^f(t) 2N . In particular, there is to > such 
that the norm is < 1/2 for t G (0, to]- Thus, Id+faEtfa is invertible for such t, and 
the previous arguments give that if K t C O t is compact then for t G (0, to], 

Ker A t n {w : suppw C K t } = {0} 

and for v supported in K t one has the stability estimate (with uniform constant C) 

\M H s f( M) < C f ||i4 t «|| ff . c -m.r-i (5 ij ) . 

We remark here that (Id +4> t Et(j)t)~ 1 can be constructed by a Neumann series, and 
thus ultimately our whole argument is completely constructive. 

In our setting we start with an ambient manifold X with equipped with a function 
x with non-degenerate level sets near the value, let x c — x + c (c near 0), let 
M c = { x c > 0}, identify a neighborhood of Y = {x = 0} with Y x (— S, 8)%, 
and have a family of operators B c £ $™ ,; (M C ) with Schwartz kernel localized in 
x < c (in both factors), where cq > is small. We further have a fixed set 
O C X with compact closure, K C O compact, and a function / continuous on 
[0,(5) with /(0) = 0, such that on O n M c , x c < f(c). In order to analyze the 
B c as c — > 0, we regard these instead as operators on M — {x > 0} by letting 
A c = B c $* c , <& c {x,y) = (x + c,y) which maps M c to M . The operators A c 

obtained by this procedure (with the parameter being c rather than t), together 
with the corresponding translates O c and K c of O n M c and K n M c satisfy all the 
requirements of the previous paragraphs, and thus conclusions apply, which, when 
translated to B c give that for sufficiently small c 

Ker B c D{w£ H^ r (M c ) : supp w C K n M c } = {0} 
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and for v G H^ r (M c ) supported in K one has the stability estimate (with uniform 
constant C) 

M\ H: ,r { M- c) < C\\B c v\\ H s- m , r - l(M - c) . 

Further, in our setting, the operators B c are in fact of the form 

B c = x- x e- F I^L c e F l^ G ^'"(Me), 

so we in fact obtain that for sufficiently small c 

KerL c n {e~ r/x ' : w G H^ r (M c ) : supp w CKf] M c } = {0} 

and for v G e F / Xa H^ r (M c ) supported in K one has the stability estimate (with 
uniform constant C) 

W e ~ F/Xcv \\H;s ( M- c) < C\\e- F /^L c v\\ m +^ (M - c y 

Notice that this is an exponentially weak estimate at dM c , i.e. at x c = 0, but the 
exponential factor is immaterial in x c > 0. Notice also that if v G H S (X), say, then 
for F > its restriction v to M c is in e F / Xc H^ r (M c ) for all r, i.e. the results are 
in fact applicable to v. 

3. Proofs 

Suppose first that X is a domain in (X,g), p G dX, and dX is geodesically 
strictly convex at p (hence near p). That is, with p a boundary defining function 
of X, we have (with G the dual metric, and metric function) that for covectors 

P e t;x \ o, 

(H GP )((3)=0^(H 2 G p)(l3)<0. 
In particular, by compactness of the unit sphere and homogeneity, there is a neigh- 
borhood Uq of p in X and Co > 0, 8 > such that for covectors /3 G X \ o, 

\(H G pW)\ < SGid) 1 ' 2 => (H gP ){P) < -C G((3). 

We then want to define a function x near p such that x(p) = 0, the region x > — c, 
p > 0, is compact for c > small, and the level sets of x are concave from the side 
of this region (i.e. the super- level sets of x). By shrinking Uq if needed, we may 
assume that it is a coordinate neighborhood of p. Concretely we let, for e > to 
be decided, an with |.| the Euclidean norm, 

x{z) = -p{z) - e\z-p\ 2 ; 

then x > —c gives p + e\z — p\ 2 < c and thus p < c; further, with p > this gives 
\z— p\ 2 < c/e. Thus, for c/e sufficiently small, the region x > — c, p > 0, is compactly 
contained in U . Further, for (3 G T* q X, H G x(j3) = -H G p(j3) - eH G \. -p\ 2 , so 
H G x = implies \H G p\ < C'eG 1 ^ 2 , so with 5 > as above there is e' > such that 
for e G (0, e'), H G x — in U implies \H G p\ < SG 1/2 , and then, for e < e', 

H 2 G x = -H 2 G p-eH 2 \.-p\ 2 > (C - C"e)G. 

Thus, there is e > such that for e G (0,e ), H%x > (C /2)G at T*X when 
H G x vanishes. Thus taking Co > sufficiently small (corresponding to eo), we have 
constructed a function x defined on a neighborhood U of p with concave level sets 
(from the side of the super- level sets) and such that for < c < Co, 

O c = {£ > -c} n {p > 0} 

has compact closure in [/ nl. 
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From now on we work with x c = x + c, which is the boundary defining function 
of the region x c > 0; we suppress the c dependence and simply write x in place of 
x c . For most of the following discussion we completely ignore the actual boundary, 
p = 0; this will only play a role at the end since ellipticity properties only hold in 
Uo and we need / to be supported in p > 0, ensuring localization, in order to obtain 
injectivity and stability estimates. Thus, completing x to a coordinate system (x, y) 
on a neighborhood U\ C U of p, for each point (x, y) we can parameterize geodesies 
through this point by the unit sphere; the relevant ones for us are 'almost tangent' 
to level sets of x, i.e. we are interested in ones with tangent vector c(Xd x + ud y ), 
c > (to say have unit length), lo e §™~ 2 , and A relatively small. 

Now, the geodesic corresponding to (z n ,u ), 7 = 7z ^ , is the projection of 
the bicharacteristic 7 emanating from {zo, gzoi^o)) = (zo,Co) (i- c - the integral 
curve of H G through this point; here we are using the metric g Zo to turn the 
vector vq into a covector) which thus satisfies (gj7)(t) = #g(7(*))j so g^(/ ° 
j)(t) = iJg/(7(i)). Thus, if / is a function on the base space X then (J^7)(0) = 
(-ffg./)(7(0),5 7 ( )(7'(0))). But Hq is homogeneous degree two in the fiber (sec- 
ond) variable of the cotangent bundle, and it is a polynomial, which shows that 
^2 (/ o 7)(0) is a quadratic polynomial in v. 

We now make this more concrete. For this, we use a fibration by level sets of a 
function x with non-vanishing differential. Letting V be a vector field orthogonal 
with respect to g to these level sets with Vx = 1, and using {x = 0} as the initial 
hypersurface, the flow of V (locally) identifies a neighborhood of {x = 0} with 
(—e,e) x x {x = 0}, with the first coordinate being exactly the function x (since 
time t flow by V changes the value of x by t). In particular, choosing coordinates 
jjj on {x = 0}, we obtain coordinates on this neighborhood such that d y . and d x 
are orthogonal, i.e. the metric is of the form f(x, y) dx 2 + h(x, y, dy), and the dual 
metric is of the form 



with /, F > 0, so (with hij denoting the metric components, so H^- is the inverse 
matrix of hij), 



F(x,y)t? + ^2H ij (x,y)rnr] j , 



dx 



= 2F{x,y)£, 



dm 

dt 



2^H l3 {x,y)r lj , 



dt _ dF 



^ 2 + H^r( x '^^ 




dt dx 



and thus 



1 d 2 x „dF , NJ . 9 „v^^F TT , . ^ 

- F(x,y)^e - Fix^^^vHVi 



which at 



dx 
dl 



0, thus ^ = 0, simplifies to 




Here we used the unit sphere for the cj-parameterization. Note that 



-^-{x,y)hik{x,y)hji{x,y)) k i 
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is positive definite by our assumptions. Thus, for geodesies we have a positive 
definite quadratic form 

a(x,y,u, 0,0) = - -g^(x,y)h ik (x,y)hji(x,y)) k iuj k LJi. 

In fact, as explained in the introduction, we mostly work in the following more 
general setting. We consider integrals along a family of C°° curves 7 x ,j/,a,u> m K™, 
(x,y, A,cj) elx M" _1 x R x §™~ 2 , depending smoothly (C°°) on the parameters, 
typically (but not necessarily) geodesies. Here R™ -1 could be replaced by an ar- 
bitrary manifold and below we make x small, so effectively we are working in a 
tubular neighborhood of a codimension one submanifold of an arbitrary manifold, 
such as X. However, since the changes in the manifold setting are essentially just 
notational, for the sake of clarity we work with R™. Further, below we work with 
neighborhoods of a compact subset {0} x K C R x x R™ _1 ; ■J x ,y,X,u(t) would only 
need to be defined for (x, y) in a fixed neighborhood U of {0} x K and for |A| < So, 
and \t\ < So, So > a fixed constant. 

The basic feature we need is that for x > and for A sufficiently small, depending 
on x, the curves stay in [0, oo) x R™ -1 . Thus, for x = only the parameter value 
A = is allowed; in our concrete setting |A| < C y/x works for suitably small 
Co > 0. However, it is convenient to use an even smaller range of A, such as 
|A| < Cqx. So concretely assume that 

lx, y ,\A°) = (x,y), l' x ,y,\, u (ty = (A,w), 
7x,j/,a,o;(*) = 2 (&(x,y, \,u,t),(3(x,y, X,u>,t)), 

and 

a(0,y,0,uj,Q) > 2C > 0, 

with a, f3 smooth. This implies that if K C R™ _1 is compact, then for a sufficiently 
small neighborhood U of {0} x K in R™ (with compact closure), and for A and t 
sufficiently small, say |A|, \t\ < Sq, where <5 > 0, one has 

a(x, y, A, oj, t) > C > 0. 

One may assume that x < So on U. Thus, writing ^(t) — (x'(t), y'(t)), 

x' = x + \t + t 2 [ (1 - s)a(x, y, A, u, s) ds > x + Xt + Ct 2 /2, 
Jo 

so if \t\ < So, (x,y) e U, |A| < <5 then 

(3-1) x>> -(*+-) +(x--). 

Thus, for |A| < V2C^/x (and |A| < <5o), x' > 0, i.e. the curves remain in the 
half-space x' > at least for \t\ < So- Further, if we fix xo > 0, then x' > x 
provided \t + ^\ > ^/2x /C and \t\ < S , thus when |A| < CoXo and |A| < So then 
x 1 > xo provided |t| > yj2xo/C , \t\ < So- Assuming x < xq and taking 

Xo sufficiently small so that ^fx + ^/IxojC < S , we thus deduce that the curve 
segments 7 :e ,; / ,a,w|(-<5o,<5o) are outside the region x' < x for t outside a (fixed!) 
compact subinterval of (— So, So)- From now on, by 7 we mean the restriction 
Jx,y.\,ui\(-8 ,s )> an d we everywhere assume that the functions we integrate along 7 
are supported in x' < x /2, so all integrals are on a fixed compact subinterval. 
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Note that in the case of geodesies, as discussed above, a is a quadratic polynomial 
in u>; this will be of use when the ellipticity of the boundary principal symbol is 
discussed. 

Before we proceed, we discuss the blowup of a space around a submanifold. Here 
we work locally on say R™, thus the submanifold can be taken to be given by w' = 0, 
where we write w — (w' , in") e R fe x M m ~ fc . Then blowing up R m_fc = {w' — 0} in 
R m amounts to introducing cylindrical coordinates around it, i.e. the factor R m_fe 
(the cylindrical 'axis', though higher dimensional) is unchanged, while on K^, one 
introduces 'polar coordinates' (\w'\, t^tt) G [0, oo) x S* -1 , thus one replaces R m by 

[R m ;R m - fc ] = [0,oo) P x S^ 1 x M"7 fc ; 
altogether one has 'coordinates' (the quotes are due to the spherical factor) 

r=\w% 0=^-, w", 
\w'\ 

with the equalities holding outside r = 0. The new boundary 

ff = {0} x S^ 1 x R m - k 

is called the front face. Further one has a blow-down map $ : [R m ; W n ~ k ] — > 
R m which is smooth, namely (r,8,w") i-> (r8,w"), but is not invertible at r — 
although it restricts to a diffeomorphism [R m ;R m - fe ] \ ff -> W n \ W n - k . We 
refer to the Appendix of [lOj for a concise but more detailed description, and for 
further references. Note that the effect of this blow up is to distinguish directions 
of approach to the submanifold being blown up, IR m ~ fc ; curves c : [0,e) — > W n 
with c(0) <E M. m ~ k and c'(0) not in the tangent space of R m ~ fc lift to (i.e. using 
the diffeomorphism property away from M. m ~ k , can be identified with) curves c 
in [K m ;K m - fc ] with c(0) € ff, and two such curves c 3 with ci(0) = c 2 (0) satisfy 
g^O) = 5 2 (0) if and only if c[(0) - c 2 (0) is tangent to R m - k . (This says that 
invariantly ff is the spherical normal bundle of W. m ~ k in M m , i.e. the quotient of its 
normal bundle minus its zero section by dilations.) 

Let X = R x x R™~ 1 , SX = XxRx §™~ 2 . In our setting, as we show momentarily, 
we start with the map 

(3.2) T+:SXx [0,oo)-> [Xx±;diag], T+(z, v,t) = j z , v {t), 

being a diffeomorphism near SX x {0}. More precisely, T + is defined on U x 
(—So, So) x §™~ 2 x (— <5o, <5o), and this map is a diffeomorphism onto it range when 
restricted to a neighborhood of SX x {0}. To see this, note that the diagonal is the 
submanifold z — z' — of X x X, so nearby one can use coordinates z — z' £ R" (the 
analogue of w' above) and z e R" (the analogue of w" above). Thus, coordinates 
on [X x X; diag] are given by z, \z' — z\ and rpEfi i an d a simple calculation shows 
that at t = 0, one has jfrEfr = (with the norms being just Euclidean norms), 
which proves that T + as in ( |3.2[ ) is a diffeomorphism near SX x {0}. Similarly, 

(3.3) r_ : SX x (-oo,0] [X x X; diag], T_(z, v, t) = 7*,„(f) 5 
is a diffeomorphism near SX x {0}. 

Remark 3.1. The analogous results would work with X in place of X. Then 

is the sphere bundle of X, i.e. TX \ o quotiented out by the R + -action. If we have 

a Ricmannian metric we could take this to be the unit sphere bundle with respect 
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to this metric, but any other choice of a transversal to the dilation orbits in the 
tangent space of X works, such as the unit sphere bundle with respect to another 
metric, or indeed (locally, in the region of interest) the space of tangent vectors of 
the form A d x + ujd y , where oj e §"~ 2 , considered above. 

We now reduce S > if necessary so that T + is a diffeomorphism on U XtV x 
(— S 0} 5 )\ x §™~ 2 x [0,(5 )t, and analogously for T_; we assume this from now on. 
(Note that in Ao we could allow an arbitrary interval with compact closure for this 
particular purpose.) 

Our inversion problem is now that assuming (If)(x, y, A, uj) = J R f(j x ,y.\,uj(t)) dt 
is known, we would like to recover / from it. (Recall our convention from above; 
the integral is really over (— S n ,S ), and f(x',y') vanishes for x' > x /2.) It is 
occasionally convenient to assume 

(3-4) 7x,j/,-A =lx,y,X,u J (t)- 

Without this symmetry assumption, we would have two curves with a given tan- 
gent line at (x,y), so having the integral of functions along both, we would have 
additional information. (In other words, we could simply drop one of these families 
to arrive at the present setting.) 

The idea is simply to average over the family, i.e. to consider for x > 

(3.5) (Af)(x,y)= [ [ (If)(x,y,X,u)x(x,X)dXdu, 

where x IS supported in |A| < \/2C^/x. One concrete choice that achieves this 

x(x,X) = x- 1 / 2 X (X/Vx), 
with x having sufficiently small support near 0; another one is 

X(x, A) = x~ 1 x(X/x), 

where now any compactly supported x works (for sufficiently small x). We remark 
that we can allow x to depend smoothly on uj and y; over compact sets such a 
behavior is necessarily uniform since there are no boundaries in these variables. 

Remark 3.2. Here we need to recall that 7 and A are only locally defined, on some 
open set O (i.e. 7 is defined for z = (x, y) e O only, and only as long as its image 
remains in O). However, as we are only interested in applying A to distributions 
supported in O, and as the ellipticity statements we show are local in nature, this 
is not a problem. For instance, for K a fixed subset of O, one may replace A by 
t\)Ai\) where ip <G C^°(0) is = 1 on a neighborhood of K, which is now globally 
well-defined, and ellipticity statements are unaffected near K. 

For any r, we can write A as 

(3.6) A = L o I, L = M 2 o n o Mi o I, 
where 

Iiu{x,y) = I I u(x,y, X,uj) dXduj, 

and 

Mm(x, y, A, uj) = x r x(X/x)u(x, y, A, uj), (M 2 f)(x, y) = x^ r f{x, y). 
Thus, n is a push-forward map, and thus is bounded on 

ff s ([0,oo) x MJ 1 - 1 xRx §™~ 2 ) ->• H s ([0,oo) x HT" 1 ) 
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for all s > 0, i.e. 'on H s ' in brief, since such a map is bounded on H s in the absence 
of boundaries, and there are continuous extension maps from H s of a half space 
to H s of the whole space. On the other hand, as \ 1S bounded, Mi is bounded on 
L 2 ({x > 0}) while its jth derivative is bounded by x~^ times a constant, so x s Mi 
is bounded as map on H s ({x > 0}) when s > integer. Thus, 

L : ff s ([0,oo) x M"- 1 x R x § n ~ 2 ) -> aT^iTQO, oo) x M™" 1 ) 

is bounded. Further, the X-ray transform, /, is itself of the form / = II o $*, where 
$* is pull-back by the map (z,i/,t) i-> J z ,v(t)i v = (A)' 4 ')) which has surjective 
differential in view of the diffcomorphism property of r± (on the relevant set; recall 
also that we are assuming that the functions we are applying / to are supported in 
U), and II is the push- forward given by integration in t. Thus, / itself is bounded 

/ : H s ([0,oo) x i"" 1 ) H s ([0,oo) x x E x S"" 2 ). 

Correspondingly, if we show A is invertible as a map between appropriate spaces 
of functions supported near x = (as discussed in the previous section), concretely 
weighted Sobolev spaces, with domain space including H s ([0, oo) x ]R™ -1 ) and range 
space including £ _s_1 -ff s ([0, oo) x IR™^ 1 ), we obtain an estimate for / in terms of 
If when / satisfies such a support condition and lies in H s ([0, oo) x R" -1 ). 



Note that the A defined by (3.5) is certainly a pscudodiffcrcntial operator in 
x > 0; moreover, its principal symbol is elliptic if \ > with \ > near (this 
uses n > 2) - while this is well-known, we check it below explicitly in the proof 
of the next proposition. Our main task is to understand the uniform behavior of 
A to x — 0. It turns out that while A itself is not a scattering pseudodifferential 
operator, its conjugates by exponential weights are: 

Proposition 3.3. Suppose \ E C^°(R). Let x(x. A) = x _1 x(A/x). The operator 
A F = x - 1 e- F / x Ae f / x is in * E T C 1 ' for F > 0. 

The main point here regarding the exponential weights is that the Schwartz 
kernel of A itself is well-behaved near compact subsets of the front face, i.e. where 
X = x -^- and Y = ^—p 1 are bounded, but is not so well-behaved as (X, Y) — > 
oo. However, the support conditions on \ insure that X — > +oo on the support 
of the Schwartz kernel of A (with a suitable estimate), and thus the exponential 
conjugation gives exponential decay of the conjugated kernel as (X, Y) — ¥ oo, giving 
the conclusion of the proposition. 

Proof. At first work in x > 0, ignoring the limit x — > 0. Then it is standard that A 
is a pseudodifferential operator (the weights are harmless then) , but it is instructive 
to prove this in a manner that extends seamlessly to the general case. 



With r± as in ( 3.2 )-( 3.3 ) , for \ an arbitrary smooth function on SX (not neces- 
sarily dependent just on x, A) the diffeomorphism property on SX x [0, 6o) allows 
one to rewrite, with \dv\ denoting a smooth measure on the transversal such as 
|dA| \du\, 

4f(*)= E I f(z'mr-\z,z'))(r-'r(\du\dt) 
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in terms of z, z' as 

' f{z')\z'-z\- n+1 b[z,^^- \z'-z\)dz', 



(3.7) 



/ z' — z \ ( z' — z \ ( z' — z \ 



\z> -z\ J ^\ \z' - z\J V ' \z> - z\. 
where a > is bounded below - it is the change of variables Jacobian factor. The 



two terms T± are in fact identical by the symmetry assumption on 7, (3.4), so 
we can ignore T_. In particular, A is a pseudodifferential operator with principal 
symbol given by the Fourier transform of 

V - z\- n+1 b(z, = V - A- n+1 {x°)(z, |^|) 

in Z = z' — z. One can insert a cutoff <f> in \Z\ with compact support, identically 1 
near (considered as an even function on M), without changing the result modulo 
rapid decay, i.e. as a principal symbol, the result is not affected. The latter can be 
computed easily as 

e- iZ <\Z\- n+1 (x<r)(z,ZM\Z\)dZ= [°° f e- itz <(x<r)(z,Z)<P(t)dtdZ 

Jo Js™- 1 

11 ' e-^<(xa)(z,Z)^(t)dtdZ=l [ j>(Z ■ ()(x<r)(z, Z) dZ; 



2 



here 4> is the Fourier transform of 4>. Fixing C — j^yi since <f> is Schwartz, if x a 
is supported away from the equatorial sphere {Z : Z ■ C = 0}, \Z ■ (\ > c|£| on 
its support for some c > 0, and then for all N, 4>{Z ■ () < CV|£| _Ar , and thus 
we conclude that the integral is Schwartz and thus gives no contribution to the 
principal symbol. Correspondingly (by using a partition of unity), it suffices to 
consider a neighborhood of the equator and assume x<r is supported here. Then 
one can write Z = (Z'^Z^) according to the orthogonal decomposition relative 
to ( = t4t, so = Z ■ similarly for Z, and dZ is of the form a(Z^)dZ^ d6, 



\z ± \ 
1 
2 



ICP 

€ §™~ 2 with a(0) = 1 since Z 1 - = (1 - \Z^ l 2 ) 1 / 2 ^ Thus, one has 
Z^C\){x^){z. l ZH + {l-\Z^ 2 ) 1 lH)a{Z\\)dddZ\\ 

^ ((ICI^ 11 \C\))a{zW) ( [ sn Jxv)(z, z"c + (1 - \zw \ 2 ) 1/2 e) do) dz\\ 



2ICI 



Since (|CI<MZ" ICD) -> <*o in distributions as |C| -> 00, this is ICI -1 J S n-2{x^)(z,0) d6 
modulo terms decaying faster as |C| — > 00; indeed, one easily sees by expanding x u 
around Z" = that this asymptotic holds modulo 0(|£| -2 ) terms. In other words, 
the principal symbol of A at (z, Q is a constant multiple of 

(3.8) Kr 1 f(x<r)(z,Z ± )dZ ± . 



In particular, if x > 0, then as long as for each (z, £), C 7^ 0, there is Z perpendicular 
to C with x non-zero at (z,Z), then A is an elliptic order —1 pseudodifferential 
operator, in accordance with the results of Stcfanov and Uhlmann [IB]. This is 
indeed the case with our choice of x, provided n > 2. 
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We now turn to the scattering behavior, i.e. as at least one of x, x' — > 0. Note 
that from (3.1 1, on the support of x, x' > x — c x 2 , for x small. We in fact show 
below that on the support of \, X is bounded below, and X — > +oo if |y| — > oo, 
and indeed X > Ci\Y\ 2 for |y| sufficiently large, C\ > 0. Here we recall from 
Section [2] that 

X = — - Y = V ~ V ' 
x 2 x 

With K denoting the Schwartz kernel of A, as 

x- 1 - (xT 1 = X —^~ =*-,= X/{1 + xX), 

XX 1 X' 

Af has Schwartz kernel 

K b (x, y, X, Y) = x^e-^-^'^Kix, y, X, Y) 

= x- 1 e- FX ^ 1+xX ^K(x,y,X,Y). 

Taking into account the polynomial bounds on K in terms of X, Y, and x' > 
x — CqX 2 implying that X is bounded below as shown later in the proof, further 
that X — ¥ +oo as \Y\ — > oo with X > Ci\Y\ 2 , exponential decay of K b as well as 
its derivatives follows easily for F > 0. Thus, the main claim is that K b is smooth 
for (X, Y) finite, non-zero, conormal to (X, Y) = 0. 

Now, onr + (suppxx [0, So)), \x — x'\ < C\y — y'\ means that locally in this region 
x, y, \y' — y\ , &E^i > \ y ~~ v \ are coordinates on [X x X; diag] - indeed, this corresponds 
to using the transversal \y' — y\ — 1 to dilations in R n = R x >- X x Kj^ in place of 
the unit sphere \{x' — x, y' — y)\ = 1, which is indeed a transversal where y' — y is 
large relative to x' — x, i.e. in our region of interest. Further, T + [x, y, A, oj, 0) is, in 
terms of these coordinates, (x, y, A, oj, 0). In general, thus, 



A(lV (x, y, \y' - y\, ^—^) )=A(x, y, \y> - y 



- x y -y 

\y'-yY\y'-y\)) "V"'*"* a " \y> - y\' \y> - y\ 



with 



/ x — x y — y 

A[x,y,0, 



\y'-y\ \y'-y\J \y'-y\ 

so (suppressing T + on the left hand side) 

\y -y\ v \y -y\ \v-v\ 

Now, in terms of the scattering coordinates, 

, / | . . x'-X xX y'-y ~ 
\y - y\ = x\Y\, — 1 = — - , — 1 = Y, 

\y-y\ \y\ \y'-y\ 

so (suppressing T + composed with the scattering blow up map on the left hand 
side) 

A X , , ~ / , , xX - \ 
(3-10) - = ^ + \Y\A(x,y,x\Y\,p-Y). 

Similarly, 

y' -y , i./ in/' i / , x' -x y' - 



+ \y' ~y\tt(x,y,\y' 



= Y 



y -y\ v \y -y\ \y -y\ 

- / xX - ' 

x\Y\n(x,y,x\Y\,—,Y 
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and 

\y> -y\ + \y> - y\ 2 f(x, y, \y> - y\, ff±- f^L\ 

v \y -y\ \y ~y\ J 

xX 



(3.11) 



\Y\+x 2 \Y\ 2 f(x,y,x\Y\,~Y 



Thus, 



(3.12) 



dtdXdio = j(x,y, —, \Y\,y) x 2 \Y\- 1 dX d\Y\ dY 
= j(x,y, \Y\,t) x 2 \Y\- n+1 dX dY 



where the density factor J is smooth and positive, and J\ x =o — 1- Also, on the 
blow-up of the scattering diagonal, {X = 0, Y = 0}, in the region \ Y\ > e\X\, thus 
on the support of x m view of ( 3.10| ), 



X 

x,y,\Y\, ^r,Y 

are valid coordinates, with |F| being the defining function of the front face of this 
blow up (i.e. of the lifted diagonal). Taking into account the x _1 in the definition 
in x, we thus deduce that K is given by 

(3.13) e-^/(i+^) x (^ + |y|A(z,y,z|n^ 

so in particular it is conormal to the front face on the blow-up of the scattering 
diagonal, of the form p^ n+1 b 1 where b is smooth up to the front face, and with- 
out the first exponential factor it, together with its derivatives (in x,y,X, Y) has 
polynomial growth estimates as (X, Y) — > oo, i.e. the derivatives satisfy bounds 
< C\(X,Y)\ N for some C, N (depending on the derivative). Decomposing K b into 
pieces supported in, say, | (X, Y) | < 2 and | (X, Y) \ > 1 by a partition of unity, we 
show in the next paragraph that the latter is Schwartz in (X, Y) due the expo- 



nential decay of the first factor of (3.13) on the support of the second factor. On 



the other hand, for the former term, supported in \(X, Y)\ < 2, calculations as in 



(3.7) in Fourier transforming this in (X, Y) show that this term of K is indeed 
the Schwartz kernel of an element of ' , with standard principal symbol being 



given by the analogue of (3.8). Here the additional information is in the behavior 
at x = 0, but given that our operator is an element of 'f^ 1,0 , the same information 
can be obtained from computing the boundary principal symbol, which we need in 
any case. 



We use (|3.11[) to express - using 



(3.14) x' = x + Xt + a{x,y, X,tu)t 2 + 0(t 3 ), y' = y + tot + 0(t 2 ), 

where the 0(t 2 ) and 0(t 3 ) terms have coefficients which are smooth in (x, y, A, cj). 
Thus, 

x' — x Xt at 2 i 3 
= x 2 = x 2 + Is? + ^2 T ( x > y > x L i > U} > t )> 
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with T a smooth function of its arguments, so 

X=^^-\Y\(l + x\Y\f{x,y,x\Y\,xX/\Y\,Y)) 



+ a(rV)|y| 2 (l + x\Y\f(x, y, x\Y\,xX/\Y\,Y)f + x\Y\ 3 T(T^), 



and thus 



, x Mr; 1 ) ^-"(r; 1 )!^! 2 

(3.15) - ± ^r 1 = \y\ +0(x) > 

where the 0(x) has smooth coefficients in terms of x, y, x\Y\, xX/\Y\, Y. Thus, for 
fi = j e (-c,c), -c\Y\ < X - a(r^ 1 )|y| 2 < c\Y\, which shows (by the positive 
dcfiniteness of a) that X — > +00 on suppx if \Y\ — > 00, and indeed, for \Y\ 
sufficiently large, X > Ci|F| 2 for some C\ > 0. 

As already explained, this proves the proposition, since now for all N' the ex- 
ponential factor in (3.131 is < C'\(X, Y)\~ N for suitable C on the support of the 
second factor, so combined with the polynomial estimates for the derivatives of 
the second and third factors, it follows that K is smooth in (x,y), with values in 
functions Schwartz in (X, Y) for (X, Y) 7^ 0, and conormal to (X, Y) = 0, which is 
exactly the characterization of the Schwartz kernel of a scattering pseudodifferential 
operator. □ 

Remark 3.4. We now explain the form these arguments would take for a different 



scaling chosen for x- By (3.14) for A = \/xfi, with /i in a compact set near (i.e. the 
first, 0(\/x) localization used above), x' = x + s/x/it + a(x, y, \fx^i, cj, t)t 2 + 0(t 3 ) 
gives 

x' <C'{x+\y-y'\ 2 ), 

indicating that \y' — y\/x 1 ^ 2 is the appropriate homogeneous variable for analysis; 
using X — , Y — lJ -^- , this amounts to a statement that the analysis is well- 

behaved on the double space of Mazzeo-Melrose [9j when the smooth structure 
is given by the boundary defining function \fx. This is a somewhat complicated 
space with a non-commutative normal operator at infinity; there's a reduced normal 
operator after a partial Fourier transform and rescaling which is a b-scattering (or 
Bessel) type pseudodifferential operator on a half-line. This is the reason for using 
our sharper cutoff, which puts us into the more amenable setting of Melrose's 
scattering calculus, as described above. 

We now compute the boundary principal symbol of Af . Indeed, this is immediate 



from (3.13) and (3.15) which show that at x — (i.e. the scattering front face) the 
Schwartz kernel of Af is 

e- FX \Y\- n+1 X ((X - a(0, y, 0, Y)\Y\ 2 )/\Y\) = K(y, X, Y). 

As described in Section[2j for each y, K(y, ., .) acts as a convolution operator, thus 
it becomes a multiplication operator upon Fourier transforming in (X, Y) , and the 
desired invertibility H^f — > 7J s s c +1 ' r amounts to the Fourier transformed kernel, 
K(y, ., .) being bounded below in absolute value by c((^,?7))~ 1 , c > (here (£,??) 
are the Fourier dual variables of (X,Y)). Thus, we need to compute the inverse 
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Fourier transform of K(y, ., .), and find F such that the desired bound holds. Note 
also that if \ depends on y and to as discussed above, we simply have 

X ((X -a(0,y,0,Y)\Y\ 2 )/\Yly,Y) 

in the above expression for the Schwartz kernel at the front face. We have thus 
shown 

Lemma 3.5. The boundary principal symbol of x~ x e~ F / x Ae F / x is the (X,Y)- 
Fourier transform of 

K(y, X, Y) = e- FX \Y\- n+l X ((X - a(0, y, 0, Y)\Y\ 2 ) /\Y\,y, t) . 

In order to find a suitable Xi we nrs t make a slightly inadmissible choice for 
an easier computation, namely we take x( s ) = e " with v to be fixed (and 
allowed to depend on y and Y), so x(.) = Cy^e - ^' ' / 2 for appropriate c > 0. 
Thus, x does not have compact support, and an approximation argument will be 
necessary. Now, in general (for arbitrary \ which has superexponential decay so its 
Fourier transform is entire), the Fourier transform in X is 

(3.16) T x K(y,^Y) = \Y f-n e -«F\Y\* e -i*s\Y\> m - iF)\Y\), 

as follows by taking into account the effect of translations, dilations and multi- 
plication by exponential weights on the Fourier transform (the last two of which 
are closely related). Here a is a function of y and Y, as above. Substituting the 
particular \ yields a non-zero multiple of 

^|y|2-n e -ar|Y| 2 e -ia||F| 2 e -^-iF) 2 |Y| 2 /2 

= y^|y|2-« e -(2«r+^ 2 -^r 2 )|v| 2 /2 e -i(a-r^|v| 2 ^ 

Now, the F-Fourier transform of |F| 2 ~ n is a homogeneous radial (i.e. SO(n — 1)- 
invariant) function of order — 1, so it is a non-zero multiple of H" 1 , with rj the 
Fourier-dual variable of Y. Notice that this uses very strongly that we have n > 2; 
n = 2 would give a delta distribution. Thus, if the Y- Fourier transform of 

e -(2aF + ^ 2 -^F 2 )|Y| 2 /2 e -i("-/^)|V| 2 

is positive, then the Fourier transform of the product, which is given by the con- 
volution (in n) of these, is also positive, and with asymptotic behavior given by 
that of \ri\~ 1 provided the Fourier transform of the Gaussian is Schwartz. Indeed, 
if one Fourier transforms \Y\ 2 ~ n ip(y, Y), where ip is Schwartz in the last variable, 
only the behavior near Y — contributes to the asymptotics as i] — > oo, and thus 
using the Taylor series of tp, one obtains the asymptotic expansion of the Fourier 
transform as r\ — > oo as a classical polyhomogcncous function (with the expansion 
in terms of j > integer). 

So assume now that a is a positive definite quadratic form in Y and take v = 
F" x a. (so same holds for v, i.e. v is a quadratic form in Y). Thus, one has ot\Y\ 2 = 
Q(Y,Y), a quadratic form in Y. Thus, writing Q^ 1 (Y,Y) for the dual quadratic 
form, and taking y(s) = e~ s ^ 2F Q( r > r )) ; we have 

X (<J) - c(F~ 1 Q(Y, y)) 1 / 2 e ~ r ~ 1< 2( CT ^' CT ^)/ 2 . 



(3.17) 
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In view of (3.16 1-(3.17|, J r xK(y,^,Y) is a non-zero multiple of 

0;|y| 2 -« e - Q / r i y i 2 e -^\n 2 e -v(ti-iF¥\Y\ 2 12 



r -l/2^ Y \2~n e -(e + F 2 )»F- 1 \Y\ 2 /2 = F -l/2 ^y^-n^ F' 1 (e+F 2 )Q(Y,Y)/2 ^ 



where the last factor is a real Gaussian since the oscillatory factor in (3.17) becomes 
identically 1. This is Schwartz in Y for F > 0, and thus the Fourier transform is a 
positive multiple of 

(det<3)" 1/2 F ( ™~ 1)/2 (£ 2 + F 2 )~ ( "~ 1)/2 e~ rcrl( ' , '' ,)/(2( £ 2+/r2)) . 

which satisfies the requirements from the previous paragraph (positive Schwartz 
function) . 

One has to be a bit careful about the joint (£, ?7)-behavior, i.e. when £ is also 
going to infinity, and where we still need lower bounds. The Fourier transform of 

\ Y \2-n e -F- 1 {£,) 2 Q{YX)/2^ with ^ = (£2 + /^l/ 2 , ig & constant mu l t iple of 



| 7? _ C |-i (e) -(— D e -^Q- i (C,C)/(2<0^) dC 

= <o _1 / i»7/(o-c / r 1 e- ro " 1(c '' c ' )/2 dc / = <e)-V(»7/(0) 



where we wrote C = C /(£)■> an d where ip is an elliptic positive classical symbol of 
order — 1, namely the convolution of with the Schwartz function e~ F ® (->-)/ 2 . 
This assures lower bounds c((£,rf)) , c > 0, i.e. elliptic lower bounds. Indeed, this 
is immediate when (£) > e j 77 j , for is equivalent to ((£, J?)) -1 m this region in 

terms of decay rates, while ip(n /(£)) is a 0th order symbol in this region. To see 
what happens when \r]\ > e{£), notice that by virtue of the classicality in fact have 
ip(rj') = \rj'\~ <p((Q/\t)'\, v' /W\), with (p smooth near in the first argument. Thus, 
we obtain 

(0-Mv/{0) = \v\- 1 m)/\v\,v'/W\), 

which is a symbol of order —1 in |£| < C 1 77 1 , and is equivalent to ((£, 

here, completing the proof of the ellipticity claim. 

Now, if \ is not a Gaussian, but rather one has a sequence Xn in C£°(R) which 
converges to the Gaussian in Schwartz functions (notice that this does not imply 
that the Fourier transform of Xn is pointwise positive for any n, which is the reason 
we need to use the Fourier transform of the Gaussian directly), then the Fourier 
transforms converge in the appropriate spaces, which suffices to conclude that for 
sufficiently large n, letting x = Xn, the Fourier transform K, i.e. the boundary 
principal symbol, still has lower bounds C*((£, f?)) -1 , C > 0, as desired. We have 
thus proved: 

Lemma 3.6. For F > there exists x G C^°(R), x > 0, x(0) = 1, such that for 
the corresponding operator x" 1 e~ ^ • ' x Ae f '• ' x the boundary symbol is elliptic; indeed, 
this holds for all x sufficiently close in Schwartz space to a specific Gaussian. 

Hence, we have 

B = x- l e- F / x Ae F ' x € *- 1 - 
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elliptic both in the sense of the standard principal symbol (in the set of interest O), 
and the scattering principal symbol, which is at x = 0, and in particular the results 
of Section [2] are applicable. Thus, elements of the kernel of A which have support 
in the compact subset K of O is finite dimensional, and further a stability estimate 
holds on a complementary subspace of this finite dimensional subspace. Further, 
with x — x c — x + c, as discussed at the beginning of this section, the arguments at 
the end of Section [2] show that for sufficiently small c, this subspace of the kernel 
of A = A c is actually trivial, and one has a stability estimate in M c = {x c > 0} for 
functions supported in K. Thus, for c small, writing the support condition as final 
subscript, 

A = xe F / x Be- F/x : e F/x H^ r (M c ) K ->■ xe F ' /x H^- r \M C ) = e F/x HI+ 1 ' r+1 (M c ) 
satisfies estimates 

Wf\\eF/*H& r (M c ) K < C W A f\\et/*m + 1 - r+1 (Mc)- 

In particular, if one is willing to give up polynomial weights as unimportant at the 
cost of losing S/x in the exponential weight, S > 0, and one uses that for s > 0, 
H S (M C ) C H°f(M c ) for r < -2±1 while for r > -5±1 + 2s, H S J{M C ) C H S {M C ), 



with continuous inclusion maps, see (2.2)- (2.3), we have 



M +s)/x H *( Mc ) K < C\\Af\\ e r/* H s+i( Mc y 



Now, using the decomposition A = L o I of A in (3.6), and the boundedness 
statements following it, we have for all F > 0, 

when / € H s+1 (M c )k ■ For the convenience of the reader, we summarize all the 
maps for F > and with r G K, S > arbitrary, in a commutative diagram: 

e F ' x H s s f{M c ) K — e F / x H s s + 1 > r+1 {M c ) — e F / x H° i r (M c ) 



H S+1 {M C ) K H S+1 {M C ) ^ F+S ^ X H S {M C ) 

with all vertical arrows inclusion maps, G the inverse of A on the range of A, and 
with 

H S +\M C ) K — L_> H°+\MoJ — ^ H S S +\M C ) 
being the lower left composite map L o I written out in detail. Indeed, note that 
even s = — 1 is allowed with the inclusions we stated; factoring the first inclusion 
via H S+1 {M C ) -> e F / x HI+ 1 - r {M c ) -> e F / x H° i r (M c ) proves it since s + 1 > 0, while 
using ( |2.4[ ) (with — s in place of s) gives the last inclusion map. In combination, 
this completes the proof of the main theorem for s > — 1 in the notation here (thus 
s > for the notation of the main theorem), with F replaced by F + S - as both 
F > and 8 > are arbitrary, this means that the original statement is proved. 
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